Disorder effect on the localization/delocalization in incommensurate potential 
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The interplay between incommensurate (IC) and random potentials is studied in a two-dimensional 
symplectic model with the focus on localization/delocalization problem. With the IC potential only, 
there appear wavefunctions localized along the IC wavevector while extended perpendicular to it. 
Once the disorder potential is introduced, these turn into two-dimensional anisotropic metallic states 
beyond the scale of the elastic mean free path, and eventually becomes localized in both directions 
at a critical strength of the disorder. Implications of these results to the experimental observation 
of the IC-induced localization is discussed. 
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The localization/delocalization of the electronic wave- 
functions is one of the most fundamental problems in 
condensed matter physics. In the presence of the dis- 
order, the interference among the scattered waves leads 
to the quantum correction to the Boltzmann transport 
theory and eventually to the Anderson localization [l|. 
This is most prominent in one-dimension, where all the 
states are localized for any finite strength of the disorder 
potential. There is another type of localization, which 
is induced by the incommensurate (IC) potential as first 
discussed by Aubry-Andre(2, yl- In this case, even in 
one-dimension, there is a critical strength of the IC po- 
tential at which the localization/delocalization transition 
occurs for all the energies simultaneously. This suggests 
that these two types of localization are different in nature 
even though the quantum interference is the key issue in 
both cases. 

It is rather easy to find the disordered system showing 
the Anderson localization, while there has been no exper- 
imental report on the IC-induced localization. Charge 
(spin) density wave system |j, |^, hehcal magnets [g, 
quasi-crystal 0, Q , and two-dimensional electrons under 
magnetic field (Hofstadter problem) Q are the represen- 
tative candidates for it. In the realistic systems, however, 
there is always some degree of disorder and it is crucial 
to see the effect of the disorder potential in addition to 
the IC potential, which we address in this letter. 

When one considers the IC potential with one wavevec- 
tor qic, which we take along x-axis, the translational 
symmetry along the two directions y, z perpendicular to 
qic remains intact, and the momenta fcy, kz are well- 
defined. Therefore the localization occurs only along x- 
direction and the wavefunctions are extended along y, z 
directions. With the disorder, one possible scenario is 



that the two-dimensional metallic state induced by IC 
potential is localized by infinitesimal disorder due to the 
Anderson localization. However, this separation into one 
and two dimensions breaks down when the disorder is 
introduced. Namely the scatterings by random potential 
reduce the spatial anisotropy. Then a keen question is 
the localization/delocalization in the presence of both IC 
and disorder potentials. Therefore, the second scenario 
is that the three-dimensional metal is recovered by the 
disorder potential. Similar problem arises for the two- 
dimensional symplectic models which remain metallic for 
weak disorder strength. 

In this paper, we consider a two-dimensional model 
with spin-orbit interaction. This model belongs to the 
symplectic universality class and shows a metal-insulator 
transition at a finite critical strength of disorder without 
the IC potential, similarly to the three-dimensional case. 
With the IC potential but without the disorder, there 
appear localized states along qic-direction at least for 
some energy region. The Hamiltonian of this model is 
given by 
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where cj, = (cJ^^^cJ^^) is the two-component creation op- 
erator, while Ct is its hermitian conjugate annihilation 
operator. The term with S in the transfer integral in 
Eq. represents the spin-orbit interaction, U is the IC 
potential, and Wr is the random disorder potential. The 
case of [/ = has been studied by Ando [ill . When the 
Fermi energy is near the band top or bottom, the non- 
perturbative part Hq reduces to the Rashba model. In 
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FIG. 1: The renormalized localization (correlation) length, 
Ac(M,W) = \aM(W)/M, for (a) parallel (a =||) and (b) 
perpendicular [a =-L) direction to qic, at the band center in 
the system with incommensurate potential of (7 = 4.0 and 

Aic = \/5. 



cumference M for several values of the disorder strength 
W. The length of a tube is taken to be iV = 4 x 10^ sites. 
We have calculated ka{M,W) with M = 4, 5, 6, 7, 8, 
10, 12, 14, 16, 20, 24, 28, 32, 40, 48, 56, and 64, and W 
is changed from 0.0 to 8.0 with the interval of 0.2. From 
the variation of the data for 1/Aq,m, we have estimated 
its error, which is at most a few percent at Aqm < 100. 
This means that the error for A|| in Fig. Q^a) is at most 
a few percent, while it is comparable to the average for 
the data Aj_ > 2 x 10^ in Fig. ^b), as suggested also 
by their scattered behavior. We present these data for 
A^ > 2 X 10^ to show only the systematic tendency. In 
Fig.^b), A_L is monotonously decreasing with increasing 
W . On the other hand. An is nonmonotonic as shown in 
Fig. nja); in the upper panel of Fig. da), the disorder 
potential W increases from bottom to the top, while it 
is reversed in the lower panel. The lowest curve in the 
upper panel of Fig. ^a) corresponds to W = Q, where 
only the IC potential is at work and the wavefunction 
is strongly localized along qic, while it is extended per- 
pendicular to it since ky is a good quantum number. ( 
The top curve in Fig. ^b) shows that X^_m reaches the 
sample size N = 4 x 10^.) 

With finite W < 2.8 , Ay shows a nonmonotonous be- 
havior as a function of M; it decreases for M < Mc{W) 
and turns to increase for M > Mc{W). This means that 
there appears a crossover length scale Mc{W), which 
is a decreasing function of W, from localized to ex- 
tended states along qic- In this range of W and M we 
have studied, it is clear that the single-parameter scal- 
ing does not hold, and it is not clear if Mc{W) diverges 
as W^ ^ Wc > or W ^ 0. However, it is convinc- 
ing that the originally localized state along qic-direction 
turns to be extended with the disorder W. A tentative 
estimate of Mc{W) from the upper panel of Fig.QJa) is 



the case of 5" <C 1, the effective mass m* and the Rashba 
coupling a are given by m* ~ l/(2ia^) and a ~ 2Sta, 
respectively. We shall take 5" = 0.5 as in Ref. [l^ and 
focus on the energy level around the band center in the 
rest of this paper. The IC wavevector qic — (27r/Aic)e2; 
is taken along the a:;-dircction with the wavelength Aic- 
The on-site potential Wr is randomly distributed in the 
range [—W/2,W/2]. In the following we shall consider 
the IC potential with U — 4.0t and Aic = \/5a. The 
energy unit (hopping parameter) t and the length unit 
(lattice constant) a are set to be 1 hereafter. 

We have calculated the localization length XaAiiW) 
of the above model in the geometry of quasi-one dimen- 
sional tube with M-site circumference in terms of the 
transfer matrix method [Il|. We consider the cases in 
which the direction of tube is parallel (a =||) or per- 
pendicular (a =_L) to the IC wavevector qic. The sys- 
tem without IC potential is also analyzed for comparison. 
Figures n^a) andQIb) show the renormalized localization 
length Aq,(M, W) = \aMiW)/AI as a function of the cir- 
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FIG. 2: Renormalized localization (correlation) length, 
A|| (M,W) = A||m(W)/M, parallel to qic as a function of 
M/^{W) where £,{W) is the localization (correlation) length 
in the thermodynamic limit. (Inset) the localization (correla- 
tion) length ^{W) as a function of W. 
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FIG. 3: Birds-eye view of renomialized localization (correlation) length, Ac,{M,W) = XaMiW)/M, for (a) parallel (a =||) and 
(b) perpendicular (a =-L) direction to qic, at the band center as a function of W and M. In each of (a) and (b),the upper and 
lower planes represents the cases without and with the IC potential, respectively. The projections of (a) and (b) to the W-A 
plane are shown in (c) and (d), respectively. 



Mc{W = 1.2) ~ 40 (7-th green curve from the bottom 
), Mc{W = 1.6) ~ 20 (9-th orange curve from the bot- 
tom), Mc{W = 2.0) ~ 10 (11-th blue curve from the 
bottom). These values are roughly consistent with the 
relation Mc{W) oc W~^, which suggests that Mc{W) 
scales with the elastic mean free path due to the disor- 
der, beyond which the incommensurability is lost and 
the system behaves as an anisotropic two-dimensional 
system. For weak disorder, the elastic mean free path 
£{W) is roughly estimated as £{W) ~ 48a{t/W)^. The 
above formula is derived from the estimations of the 
mean free time r as t"-'^ ~ Tr{w'^)Np, the density of 
states per unite cell Np ~ l/(27ri), the averaged dis- 
order strength {w'^) ~ W'^ /12 for the rectangular distri- 
bution, and the Fermi velocity vp ~ 2to around the band 
center. This estimation semi-quantitatively supports the 
identification, MdW) ~ ((W), since i{W = 1.2) ~ 33, 
£{W = 1.6) - 19 and i{W = 2.0) - 12, for this regime 
of disorder. 



In the lower panel of Fig. ^a), on the other hand, 
W increases from the top {W — 3.0) to the bottom 
(W = 8.0). AtW ^ 3.0, MciW) is almost the lattice con- 
stant and A II shows the monotonous behavior as a func- 
tion of M for Vl^ > 3.0. This is similar to A_l in Fig.mb), 
and is the canonical behavior of the metal- insulator tran- 
sition. Therefore we show in Fig. |21 its single-parameter 
scaling analysis |12!| from the lower panel of Fig. Q^a). 
The degeneracy of the data indicating the scaling rela- 
tion 
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is clearly obtained and the resultant localiza- 
tion/correlation length £,(W) is plotted in the inset 
of Fig. El (In the upper branch, only the data for M ^ 8 
are used.) 

The divergence oi£,{W) toward Wc = 4.3±0.1 does not 
contradict with the relation ^{W) ^ \W — Wc\~'^ with the 



known critical exponent v = 2.7 ± 0.1, although we did 
not try the accurate estimate of v. Therefore we believe 
that the present metal-insulator transition belongs to the 
universality class of two-dimensional symplectic one |lCt 

Ellli. 

In order to obtain a global picture, we plot in Fig- 
ures|3Ja) andlJ^b) the birds-eye view of the renormalized 
localization length ka{M,W') = XaM{W)/M as a func- 
tion of the circumference M and the disorder strength 
W for (a) a =|| and (b) a =_L. In Figs. Etc) andEtd), 
we also put the projections of Figs.EJa) and^Ib) to the 
W-h-a plane, respectively. In each panel, two curved sur- 
faces are shown. The upper (lower) surfaces are for the 
case without (with) IC potential. It should be noted that 
the color contour represents the value of M in logarithmic 
scale, different from the colors in Figs.Handl^l in order to 
make the connection to the projected figures on the right- 
hand side. From the A|| with U = 4.0 in Figs.Eta) and 
Ofc), it is clearly seen that the state is localized, i.e., the 
red part is lower than the blue one, at both small W = Q 
and large VF > IVc = 4.3±0.1. The degeneracy occurs at 
Wc indicates the metal-insulator transition as analyzed 
in Figs.^andl^l while one can not see any degeneracy for 
smaller W . Therefore, the reversal of the M-dependence 
around W = 1.5 is attributed to the crossover at Mc{W) 
as discussed above. The behavior of Aj_ in Figs.OIb) and 
|3Jd), on the other hand, is similar to the case oiU — Q. 
Note that the critical Wc is common for both A|| and A^, 
which is reduced from Wc = 5.8 ± 0.1 to Wc = 4.3 ± 0.1 
by the IC potential U = 4.0. 

Now we discuss the relevance of these results to the 
experimental observation of the localization due to IC 
potential. Except the one-dimensionally confined sys- 
tems such as quantum wire, the electronic wavefunctions 
are extended in the perpendicular directions to the IC 
wavevector. This is almost always the case as in charge 
(spin) density wave, and helical magnets. In this case, we 
expect that the elastic scattering by impurities washes 
out the IC potential in both two and three dimensions 
beyond the length scale of the elastic mean free path 
1{W). First let us ask the following question: " Can 
one observe the temperature dependence of the resistiv- 
ities characteristic to the IC-induced localization ? " To 
answer to this question, one should remember that the 
temperature dependence of the resistivity is translated 
into that of the inelastic mean free path Lin{T), which 
is replaced by the sample size L. The quantum correc- 
tion due to the interference occurs for the length scale 



L > £{W). Namely, £{W) is the shortest length scale for 
the localization problem, and the Boltzmann transport 
theory applies for L < £{W). Therefore, if the crossover 
length scale Mc{W) is identical to i{W), as suggested by 
the analysis above, there is no chance to observe the IC- 
induced localization from the temperature dependence. 
Most probably this is the reason why there are no exper- 
imental observation of the IC-induced localization, i.e., 
it is very fragile against the disorder and is replaced by 
the usual Anderson localization. However, the incom- 
mensurability gives a strong anisotropy of the metallic 
conduction between parallel and perpendicular directions 
to the wavevector qic within the Boltzmann transport 
regime, which is the remnant of the IC-induced local- 
ization. Therefore, for example, one should look for the 
dramatic change of the Resistivity anisotropy associated 
with the commensurate-incommensurate transition. 
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